Quasi-spherical gravitational collapse and the role of 
initial data, anisotropy and inhomogeneity 



In this paper, the role of anisotropy and inhomogeneity has been studied in quasi-spherical gravi- 
tational collapse. Also the role of initial data has been investigated in characterizing the final state 
of collapse. Finally, a linear transformation on the initial data set has been presented and its impact 
has been discussed. 

PACS numbers: 



Gravitational collapse is an important and challenging issue in Einstein gravity for the 
last two decades or more, particularly after the formation of famous singularity theorems [1] 
and Cosmic Censorship Conjecture (CCC) [2]. Moreover, the final outcome of gravitational 
collapse [3] in the background of general relativity is interesting both from the point of 
view of black hole physics as well as its astrophysical implications. In fact, an extensive 
study of gravitational collapse [4 - 9] has been carried out of Tolman-Bondi-Lemaitre (TBL) 
spherically symmetric space-times containing irrotational dust to support or disprove the 
CCC. A general conclusion from these studies is that a central curvature singularity forms 
but its local or global visibility depends on the initial data. Recently, mena et al [10] 
have introduced a linear transformation on the initial data set which keeps p and a to be 
invariant at any time instant except for the initial ones and have discussed the consequences 
for the final fate of collapse. 

On the otherhand, there is very little progress in studying non-spherical collapse due to 
the ambiguity of horizon formation and the influence of gravitational radiation. However, in 
recent past, an extensive study of irrotational dust collapse has been done in quasi-spherical 
Szekeres space-time [11-13]. In this paper, we have extended the work of Mena et al 
to quasi-spherical Szekeres model and have investigated the effect of inhomogeneity and 
anisotropy on the final nature of singularity. The paper is organized as follows: a brief 
outline of the quasi-spherical dust collapse has been presented in section II. Section III deals 
with properties of inhomogeneity and anisotropy both initially and at any time instant and 
then linear transformation on the initial data set have been discussed. The paper ends with 
conclusion in section IV. 



The inhomogeneous quasi-spherical dust collapse in n-dimension is represented by Szek- 
eres' space-time with line element 
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I. INTRODUCTION 



II. QUASI-SPHERICAL DUST COLLAPSE 
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where the metric coefficients a and (3 are functions of all the n space-time co-ordinates and 
have explicit form [12,14] 



The evolution equation for R is 



e a = R l + RV ' , eP = R(t, r) e v ^-^) 



R 2 = f(r) 



F(r) 
R n ~ 3 



and v has the explicit form 



n-2 n-2 



e-»=A(r)J2x 2 + Y, B ^ +C ^ 



i=l 
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(3) 
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where F(r)(> 0) and f(r)(> —1) are arbitrary functions of r and A(r) 7 i?i(r)'s and C(r) 
are arbitrary functions of r alone with the restriction 



(5) 



i=i 



(dot and dash stands for partial differentiation with respect to t and r respectively). 
The matter density and the anisotropy scalar are given by [12] 



and 



p(t,r,xi, ...,x n - 2 ) = 



a(t,r,X!, ...,x n - 2 ) 



(n-2) 



F' + (n - l)iV 



R n ~ 2 (R' + Rv') 



n-2 



2(n- 1) 

Now integrating the evolution equation (3) we have 



RR f — RR r 
R(R' + Rv') 



(6) 



(7) 



t-U 



[n-l)VF 



A fr n -\ n-i A fR 

r— 2 F 1 [-,a,a+l,-^ r -]-R— 2 *i[-, a, a + 1, -i—^ 



n—3 



(8) 



where R = r at the initial epoch t — U and 2 F\ is the usual hypergeometric function with 

2 ^ n-3 ■ 

If i = t s (r) stands for time of collapse of the r-th shell i.e., R(t s (r), r) — then we have 



2 n-i 1 f r n-3 

ts(r) ~ U= (n-l)^ ^ 2Fl[ 2 ' °' a + 



(9) 



Further, at the time of formation of trapped surfaces (i.e., t — t a h(r)) we have R n 3 — 
F(r), so we get 



t a h(r) - U 



OrV I f r n-3 OF^^ 1 

jFil-,0,0 + 1,-^]-— r2 Fi[-,o,o+l ) -/] (10) 



{n-l)VF L 2' 



n- 1 
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The time of formation of central singularity is given by 

n— 1 n—3 

t Q =t s (0)=U+ Um (S}7F 2F 1 [la,a+l,-f^] (n) 
r -> 

For finite value of the above limit, we assume F(r) and f(r) to be in the following poly- 
nomial form near the central singularity (r = 0) 

F(r) = For™- 1 + F x r n + F 2 r n+1 + , 

(12) 

f(r) = for 2 + fir 3 + / 2 r 4 + 

Then equation (6) suggests that initial density profile pi (= p(ti,r)) is smooth at the 
centre and we have 

Pt(r) = po + pir + p 2 r 2 + (13) 

Now using the initial condition R{ti,r) — r the initial density and shear have the expres- 
sions 

_ (n-2) F' + {n-\)Fi/ 
Pl 2 r"" 2 (l+n/) 1 j 

and 



/ n-2 [{rF'-(n-l)F} + r"- 3 (rf -2/)] 
Y8(n-1) r ¥ (F + /r«- 3 ) 1/2 (l + ri/) 

Using the series expansions (12) and (13) in equation (14) and comparing terms of equal 
powers in r we have 

P J -= ( " + J "~ 2 1)( "~ 2) ^ i = 0,l,2,.. (16) 

Further, the above series expansions result (after simplification) an explicit expression for 
the time difference between the formation of trapped surface and the occurrence of central 
singularity as (upto leading order in r) 



2 -J- 

tah-ta = -F n 3 r—3 



1 u (n-l)VF 



^ 2Fl[ 2' a ' a+1 '-^ ] 



(n-l)/ //i iq\ 3 ^ /o, 

+77; ^7-7-77 2Fi[-,a + l,a + 2,-— ] 



(3n-7)F V/o F o y iL 2' ' ' F 



(17) 



The necessary condition for existence of naked singularity is characterized by t a h{r) > to- 
Hence the nature is completely specified by the initial data set X = {F, /}. It is to be noted 
that if the coefficient of r in the expression (17) vanishes then the end state of collapse also 
depends on the dimension of the space-time. 



III. A STUDY ON ANISOTROPY AND INHOMOGENEITY 



In this section we shall study some characteristic of anisotropy and inhomogeneity of the 
space-time, we shall state them in the form of proposition. 
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A. Some Results 

Proposition I: Ifa(U,r) = = &(t i} r) then the initial data set is homogeneous and the 
initial space-time is homogeneous and isotropic. 

Proof : The expression for &(U, r) is 



&(U,r,Xi, ...,x„- 2 ) = 



n-2 (n-3)[(n- l)F-rF'} 



2(n-l) 2r(l+rz/) 



o-i 



F(b-n+2rv') 



+ /(l + rv') + r (/' + ^) 



Hence a = = <7j implies 



2r(l + n/')V^ + / 
(n - 1)F - r,F' = 



i.e., 



Also <7,: = results 



F = F r 



n-l 



So using (19) we have 

r/'-2/ = i.e., / = / r 2 . 
Thus the initial data is homogeneous. 
The expression for initial density contrast is 



dp 

dt 



(n-2) 



F" (n - 2)F' 

r n - 2 (l + r^') r^-^l + r^') 2 



(18) 

(19) 
(20) 

(21) 

(22) 



+ 



{-(n - l)(n - 2)iV + r (rf - (n - 1)F) ((n - l)t/' 2 - v") } 



(23) 



which vanishes identically for the choice of F in equation (20). It is to be noted that this 
choice of F simplifies the initial density to 

(n-l)(rc-2) 
Pi g 

which is homogeneous. Hence the initial space-time is also homogeneous and isotropic. 

Proposition II: If a(t, r, x\, x n -2) = then (i) the initial data set is homogeneous, 
(ii) the space-time is homogeneous and isotropic, (Hi) the end state of collapse will always 
be a black hole. 
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Proof : From equation (7), a = implies 

~R' ~ R 

which shows that R = g(r)fi(t). But at t = U, R = r, hence we have 

R{t,r)=r»(t) (24) 
Now substituting in the evolution equation for R (i.e., eq.(3)) we have 

A' - 77^3 + fo (25) 
with F = For" -1 , / = /o?' 2 i-e., the initial data set is homogeneous. 
So the expression for density simplifies to 

_ (n-l)(n-2)F 
P[ ) ~ 2/x"" 1 

Hence the space-time is homogeneous and isotropic (by condition). 

Further, using this expression for R (i.e., eq.(24)) the time difference between the forma- 
tion of trapped surfaces (t a h) and the occurrence of central singularity (to) becomes 

t a h-t = -j^—^F^r^ 2 F 1 [i,a, a +l,-/ r 2 ] (26) 

So always we have t ah < t and the end state of collapse will be a black hole. It is to be 
noted that similar results holds also for spherical collapse [14]. 

Proposition III: If p' i = then the initial data set is homogeneous as well as the initial 
space-time is homogeneous and isotropic. 

Proof : The expression for p ■ can be simplified to 

= 2 ^ n ( r+U 2 [(1 + rvl){rF " - {n 2)F ' } 



+{rF' - (n - l)F}{(n - 2) + (n - \)rv n - rv"}\ (27) 

which shows that for p[ = we can choose F = For" -1 . Also for regularity of initial 
shear at r = demands / = for 2 . Hence the initial data set is homogeneous. For this 
choice of F and /, pi becomes constant and <r, = 0. Therefore, the initial space-time is 
homogeneous and isotropic. 

Proposition IV: If the density contrast vanishes at any time instant (i.e., 
p'(t, r, xi, £„_ 2 ) = 0), starting from regular initial data, then the space-time will 
be homogeneous and isotropic and final state of collapse will be a black hole. 

Proof : Using the condition for regular initial data (i.e., F = F r" -1 , / = /or 2 ) the 
expression for energy density simplifies to 

, (n-l)(n-2) r"- 2 (l + ri/)F) 



G 



which shows that if wc choose R = r[i(t) then p will be independent of r i.e., 
r, Xi, x n -2) = 0. It is to be noted that p will also become independent of other 
space co-ordinates. Also the choice R = rp(t) results a = 0. Hence the space-time is 
homogeneous and isotropic. In proposition II, we have shown that a = results to a black 
hole as the final state of collapse. Hence the proof. 

From the above propositions it is clear that if the space-time is homogeneous and 
isotropic then the final outcome of collapsing process will be definite (a black hole) but no 
such conclusion is possible if the space-time is inhomogeneous and anisotropic. 



B. Transformation of the initial data set 

In the initial data set X = {F, /} if we make a linear transformation [10] 

C T : 1 -> {a^F, af} (a > 0, a constant) (29) 

then the scale factor R changes as R — > a x ^ 2 R while (p, a) are invariant under this 
transformation but not the initial density and shear i.e., (pj,<7j). In other words, there are 
classes of values of F and / for which we have the same values of p and a but the initial ones 
(namely pi and <7j) are different for each set of (F, /). So it may happen that for one initial 
data set (F, /) we have naked singularity as the final end state but choosing 'a' appropriately 
we have the transformed (F, /) for which we may black hole as the final outcome of the 
collapsing process but p, a remain same for both the initial data set. This has been shown 
numerically for spherical collapse by Mena et al [10]. So we may conclude that energy 
density and anisotropy scalar at any time instant can not uniquely characterize the nature 
of the final state of collapse. But for the initial nature of the space-time no such conclusion 
can be made. In fact for a given initial inhomogeneity and anisotropy the initial data set can 
be specified uniquely and hence final outcome of the collapsing process may be characterized. 



IV. CONCLUSIONS 

In this paper, we have studied properties of anisotropy and inhomogeneity of the 
space-time both initially and at any epoch in the context of gravitational collapse. We have 
obtained two distinct features from these studies namely: (a) the matter density and shear 
scalar at any instant remain invariant under the linear transformation (29) and hence can 
not characterize the nature of the final singularity, (b) the initial energy density and shear 
scalar are not invariant under the linear transformation (29) but together can uniquely spec- 
ify the initial data set and hence they may characterize the final outcome of the collapsing 
process. This is a very peculiar feature of the physical process because physical parameters 
like energy density, shear must have consistent behaviour for all time, they should not have 
different characters initially and at any epoch. It seems that the trouble is with the ini- 
tial choice of the scale factor. We shall discuss the role of initial choice in a subsequent paper. 
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